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Abstract General relativity eonsiders Dixon's theory as the standard the- 
ory to deal with the motion of extended bodies in a given gravitational back- 
qh| ground. We discuss here the features of the "reaction" of an extended body 

^ . to the passage of a weak gravitational wave. We find that the body acquires 

a dipolar moment induced by its quadrupole structure. Furthermore, we de- 
rive the "world function" for the weak field limit of a gravitational wave 
background and use it to estimate the deviation between geodesies and the 
, world Hues of structured bodies. Measuring such deviations, due to the exis- 

• fence of cumulative effects, should be favorite with respect to measuring the 

' amplitude of the gravitational wave itself. 
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1 Introduction 

The equations of motion for an extended body in a given gravitational back- 
ground were deduced by Dixon l,2,3,4,5j (hereafter the "relativistic ex- 
tended body model," or simply "Dixon's model") in multipole approximation 
to any order. In the quadrupole approximation they read 

-^R^uc.l3U''S"f^ - ^J°''^''^Ra0jS''" = i^f'^P'")^ + F^'i"'^^)'^ , (1) 
2p[A'[/H + 1 J^^Tl^i?'']^^^ , (2) 

where = mU^ (with Up - Up = —1) is the total four-momentum of the 
particle, and 5^'' is a (antisymmetric) spin tensor; U is the timelike unit 
tangent vector of the "center of mass line" Cjj used to make the multipole 
reduction, parametrized by the proper time tij. The tensor J^/'t* is the 
quadrupole moment of the stress-energy tensor of the body, and has the 
same algebraic symmetries as the Riemann tensor. In this paper we limit our 
considerations to Dixon's model under the further simplifying assumption 
[Sl[7] that the only contribution to the complete quadrupole moment J°'Pi^ 
stems from the mass quadrupole moment Q"'', i.e. we write 

= -3i7^°Q^l[TC/p^l , Q'^^Upp = . (3) 

In order the model to be mathematically consistent the following addi- 
tional condition should be imposed [I] to the spin tensor 

S^^'Up, = , (4) 

to ensure the correct definition of the various multipolar terms. 

It is also convenient to introduce the spin vector by spatial (with respect 
to Up) duality 

~ 2'7a^'' Up S-/S , (5) 

where rjap-ys — ^Z—g^a/i-yS is the unit volume 4- form and ea/3-ys (eoi23 = 1) is 
the Levi-Civita alternating symbol, as well as the scalar invariant 

= ^5^.5^'' = S^S^^ , (6) 

which is in general not constant along the path. 

There are no evolution equations for the quadrupole as well as higher 
multipoles as a consequence of the Dixon's construction, so their evolution is 
completely free, depending only on the considered body. Therefore the system 
of equations is not self-consistent, and one must assume that all unspecified 
quantities are known as intrinsic properties of the matter under consideration. 
Moreover, the assumption that the considered body is a test body means 
that its mass, its spin as well as its quadrupole moments must all be small 
enough not to contribute significantly to the background metric. Otherwise, 
backreaction must be taken into account. 
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We investigate here how a smah extended body at rest reacts to the 
passage of a single plane gravitational wave (GPW) in the weak field limit. 
The body is spinning and also endowed with a quadrupolar structure. The 
latter property allows us to consider this analysis as a good model for a 
bar antenna. The present work thus extends previous results by Mohseni 
and Sepangi [8] concerning the motion of purely spinning particles (i.e. with 
negligible quadrupolar structure) in a GPW background. 

2 Motion of extended bodies in GPW spacetimes 

Consider the metric of a single gravitational plane wave propagating along 
the X direction of a coordinate frame, with "+, x" polarization states, written 
in the form 

ds^ = -dt^ + dx^ + (1 - h+)dy'^ + (1 + h+)Az^ - 2hy,dxAz , (7) 

where /i+,x = h^^^it — x) are functions only oi t — x because of vacuum 
linearized Einstein equations. Let the metric functions have the following 
form 

hj^ — sin Lo{t ^ x) , — cos Lu{t — x) . (8) 

A wave with linear polarization is characterized by = or Ax = 0, 
whereas the condition A^ = ±Ax corresponds to a circularly polarized wave. 
It is also useful to introduce the "polarization angle" of the wave tan^/; = 
Ax/A^. Throughout the papers all tensorial quantities are assumed to be 
decomposed along the coordinate frame {bq} = {da}- 
The geodesic of this metric are given by [9] 

t^(gco) = ^[(m' + / + E')dt + (m' + / - E')d.] 

+ + h+) + f3hx]dy + - h+) + ahx]d, , (9) 

where a, [3 and E are conserved Killing quantities, — 1, 0, —1 correspond- 
ing to timelike, null and spacclike geodesies respectively, and 

/ = a^{l + h+) + (3^{1 - h+) + 2a(3hx . (10) 

The corresponding parametric equations of the geodesic orbits arc then easily 
obtained 

t(A) =EX + to+ x{X} - xo , 

x{X) = ifi^ + a^+p^-E^)^- ^[(«' - P')A+ coscu{EX + to - xo) 

—2af3Ax smuj{EX + to — xq)] + xq , 
y{X) =aX + yo 

Xo) - PAx sinw(^;A + to- xq)] , 



Xo) + aAx smui{EX + to - xo)] , (11) 



— \aA+ cos uj (EX + to — 

LUE 

z{X) = pX + zo 

H — ^ [l3A+ cosLj{EX + to- 
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where A is an affine parameter and Xg are integration constants. 

For later convenience we also evaluate the world function, which is defined 
for any given background metric as half the square of the spacetime distance 
between two generic points xa and xb connected by a geodesic path [TU] 

n{xA,XB) = -j^ g^.[x«(A)]^^dA^ r2^B , (12) 

where a:"(A) satisfy the geodesic equations and the afhne parameter is such 
that a;"(0) = x^ and x°' {!) ^ x% . We obtain 

+^[{yB - Va)^ - {zb - za)'] ^ 

2u tB ~ Xb - [tA- xa) 

Ay, smuj{tB ~ Xb) - smuj{tA - Xa) 
iVB - yA){zB - za) 7- T , (13) 

UJ tB- Xb - [tA - XA ) 

where 

^?AB = \'nc.(i{xA - Xb)"(x^ - Xb)P (14) 

is the well known expression of the world function in the Minkowski space- 
time. To our knowledge the expression p3p for the world function cannot be 
found in the literature. 

It is in general a hard task to solve the whole set of equations (dJ-Q 
in its complete generality, without imposing convenient restrictions on the 
form of the various quantities involved, e.g. constant frame components of 
the spin and quadrupole tensors, center of mass line moving along a given 
orbit, total four- momentum of the body aligned with a given direction [llj . 
However, since we are interested in a weak field solution it is sufficient to 
search for changes in the mass parameter m of the body, its velocity and 
total four-momentum as well as spin and quadrupole tensors which are of 
the same order as the the metric functions h^^x- 

m^nio + m , U ^ Uq + U , P = Pq + P , 

^M,. ^ 5-^"^ + 5-^^ ^ QM^ = Q^-^ + QM^ ^ (15) 

where the quantities with the subscript "0" are the flat background ones, 
i.e. those corresponding to the initial values, before the passage of the wave. 
By taking advantage from that we are able to solve the problem in its most 
general form. The background quantities are solutions of the zeroth-order set 
of equations 

^ - , = 2p[^[/o"l . (16) 



druo druo 

Assume that the body is initially at rest at the origin of the coordinates, i.e. 
the associated world line has parametric equations 



t = Tuo, x{tuo)=0, y{TUo)=0, z{tu„) = 0, 



(17) 



where tu^ denotes the proper time parameter, and hence the unit tangent 
vector reduces to C/q = dt- The set of equations (fTHll is fulfiUed for instance 
by taking Pg = 'm-oUo, implying that the components of the background 
spin tensor S'q" remain constant afong the path. The conditions Q simply 
give 5*0° = 0. The remaining components are related to the background spin 
vector components by 

S12 q3 cfl3 q2 023 q1 /-I Q\ 

Q — Dq , »-'0 — '-'0 ' '-'0 — '-'0 • \^°) 

Note that the quadrupolar quantities do not enter the flat spacetime equa- 
tions, since they are coupled to the Riemann tensor. Therefore, we can as- 
sume the components of the background quadrupole tensor Qq" to be also 
constant with respect to the coordinate frame, without loss in generality. The 
conditions dS]) give Qg" = 0. 

The conditions ensuring that the body does not perturb the metric can 
be formulated by using the natural length scales associated with it, i.e. the 
"bare" mass toq, jS'g |/too and HQo^l/mo)^/'^ . 

Let us turn to the general system of equations ([I])-(l4]). The first order 
set of equations is obtained after substituting in it the perturbed quantities 
(|15p . by retaining terms up to the first order in the metric functions. Let 
the "center of mass line" be generic after the passage of the wave, i.e. with 
timelike unit tangent vector U given by 

[/ = 7(9t+^"9„) , 7= (l-i>2)-i/2 ^ i}^ = 5abi>''i>'' , a, 6= 1,2,3. (19) 

According to the decomposition (|15p the linear velocities along the spatial 
axes are first order in /i+.x, so that 7 ~ 1. Analogously, let the total four- 
momentum P = mUp be generic, with 

Up = 7p(ft + i>;da) , 7p - (1 - i^pr'^^ , ^p = Sabi>;^'p ; (20) 
also in this case 7^ ~ 1. The supplementary conditions ^ imply that 

S^^ = Sq^p — S^Vp , 5^"^ = —Sf^Vp + Sq^p , S'^^ = S^Vp — S^Vp . (21) 
The spin vector is orthogonal to Up from its definition ([5|) , with components 

5" = Slv^ + Slv^ + Sp^ , = Si + ^ 

52 ^ Si - , 53 ^ ^3 ^ ^12 (22) 

Finally, the conditions ([31)2 on the components of the first order quadrupole 
tensor give 

qOO = , = -Ql'^vl ~ Ql'^vl ~ Ql'^vl . (23) 

It is useful to introduce the rescaled dimensionless quantities (for both 
zeroth order and first order terms) 
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Note that the components of the zeroth order quadrupole tensor all have 
definite sign, whereas the components of the spin vector are allowed to have 
both signs, corresponding to a body spinning up/down. 

Consider first the equations of motion ([1]) . The spin force and the quadrupole 
force as defined in Eq. ([1]) turn out to act on the orthogonal two-planes y—z 
and t — X respectively: 

^(quad) ^ [^^^ ^^^j^ _ 2 tan V^/i+q23] (Q, _ a,) , (25) 

where we have introduced the quantity /o = q^-^ — q^^ . Eqs. (P) thus reduce 
to the following set 

_ ^ = ^(quad)O ^ TOo^ = i^(quad)l ^ 

dr(7 ' dru 

mo-^ = F(«P-)2 , mo-r^ - F(^p-)3 ^ (26) 
dru dru 

whose solution is straightforward: 

m = [/o^+ sincjTc/ + 2q'^^Ax {cosljtu - 1)] , 

vl ^ — - 4r'?o^ + "^i ^ ^7 [/o^+ sinwT;/ + 2ql^Ay, cosuitu] + ci , 
= — - [ctq^x svaujTu + crQA+ coswr^/] + C2 , 



f 2 
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[ctq^x sinojTt/ — (jIA^ cos wt[/] + C3 , (27) 



where the initial condition m(0) = has been imposed, and ci, C2 and C3 
are integration constants. Eq. ([27|)i implies that the zeroth order quadrupole 
tensor is responsible for the mass of the body to vary after the passage of 
the wave. 

Consider then the evolution equations ^ for the spin tensor. By using the 
supplementary conditions (|2ip they give three algebraic relation between the 
spatial linear velocities of U and of Up plus three evolution equations 
for the spatial components 5" of the spin tensor to be integrated together 
with the initial conditions (t"(0) = 0. The corresponding solution turns out 
to be given by 

= i[2((fT2)2_(ai^)2) + /o]^^sincjrt, + i(2a2aij+923)^^(coswrj/-l) , 
+ ql^) + Ay,{al ~ alqf )] sinwr^/ 
-i [2A^{<jI<jI + ql^) + A+<jI{2 + /o)] (cosc^rc; - 1) , 
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= ^/o^x s'mujTu + qQ^A+{cosujTu - 
= i(Axcro - ql^A+) (cos ujTu - 1) - ■ 
= i(Ax(7^ - q^'^A+) (cos ujTu - 1) - ■ 



/o)] (C0SCJT[/ - 1) , 
1), 

^{A+(tI ~ ql'^Ax)sinuJTu , 

i(A+f73_ql3A^)sinCJTt; , (28) 



where the constant ci , C2 and C3 have been conveniently fixed so that i>°' (0) = 
0, implying that 

i>l = [foA+ sinLJTu + 2ql^ A ^ (cos ujTu - 1)] - (4<J^ + qo^)A^ , 
i>p = [cr^Ay^ sinujTu + alA+(cosuJTu - 1)] 
+ 1 [2^x (aoVg + + A+a'jo] , 

l>p = i [ctqAx sincJTjy - crQ^+ (cos WT(7 - 1)] 

+ ^[2A^(alal+ql^)-A+alfo] • (29) 

Eqs. ([281)4,5,6 thus imply that even if initially absent the spinning structure 
will be acquired by the body during the evolution due to its quadrupolar 
structure. 

Finally, the nonvanishing components of the first order quadrupole tensor 
can be easily obtained from Eqs. (I23p . whereas (7° is given by Eq. ([22]) . 

The modification to the initial trajectory (|17|) of the body after the 
passage of the wave is easily obtained by integrating Eqs. ([28]) 19^. since 
i>°' — da;°/dr[/, together with the initial conditions x'^(O) = 0: 

t = Tu , 

Lox = -i[2((a2)2 - (alf) + h]A+(cosujTu - 1) 

+ ^[2o-oc^o + Qx)^]A>c (sin LJTy - ujtu) , 
^y^l [^+(^o'^o+'?o')+^x(^o-^o'Zo')] (C0SUJTU~1) 

-i [2Ax (dldl + ql>^) + A+al(2 + fo)] (sinojTu - ojtu) , 

'^^ = [^+('^o'^o +90^) + ^x(^To - c^o^o^)] {COSUJTU - 1) 

-i [2Ay,(alal + q^) + A+al(2 - /„)] (sinwrcj - ujtu) . (30) 
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3 Discussion and general results 

In the expressions (|30p we can distinguish a harmonic part and a secular one; 
the latter is given by 

w j/scc = \ [2Ax {alal + Qo") + ^+^^^(2 + /o)] ^^tc/ , 

tJZsec = ^ [2Ax(fT,icrg +90^^) +.4+(72(2- /o)] tJTc/ , (31) 

and it is important in view of "cumulative" effects which it originates. 

It is now interesting to evaluate the deviations between the geodesic world 
line of a reference test particle (always) at rest at the origin of the coordinates 
and that of the structured body under consideration here, initially (only) at 
rest at the origin of the coordinates, i.e. with unit tangent vector initially 
aligned with the geodesic one. The deviations can be obtained by evaluating 
the world function (|13|) at the same coordinate time tu^ — t — tjj, where 
x'^ = x"\rug and x% = x"\ri, are given by Eqs. pT|) and respectively. 
The world function is given by 

f^AB = ^SabX'^x'' , (32) 

by neglecting higher order terms in the amplitudes A+ and j4x of the gravi- 
tational wave. Apart from oscillating contributions, the dominant one is 

SIab ^ ^ [a^scc + Vscc + 2soc] > (33) 
so that the deviation turns out to be 

dAB = ^/2\nAB\ ^ V X^^e + Vscc + ^soc • (34) 

Moreover, assuming ~ Ay ~ h, from Eqs. ([51]) it follows that the devi- 
ation dAB hTjj, the constant of proportionality depending on the proper 
structure of the body. It is therefore clear that cumulative effects arise and 
can be used to get better conditions to measure (indirectly) the amplitude 
of a gravitational wave. 

The complementary cases of an initially non-spinning body with a given 
quadrupole structure as well as that of a purely spinning body with a negli- 
gible quadrupole structure will be considered in detail in the Appendix. The 
results concerning two simple cases are listed below. 

Consider first the case of an initially axially symmetric (about the z-axis) 
non-spinning body described by the zeroth order quadrupole tensor 

q^" = diag [/o/3, /o/3, -2/0/3] , (35) 

where the tracefree property has been assumed in analogy with the classical 
(non-relativistic) case. Eqs. (|4ip imply no secular deviation, since 

c/ab = ^|cosu;tc/-1| . (36) 



9 



Note that for a non-spinning body secular deviations can arise only for non- 
vanishing nondiagonal components of the quadrupole tensor, as from Eqs. 

(ED). 

Consider then the case of a purely spinning body along the y-axis for in- 
stance, with vanishing quadrupole tensor. Eqs. ()45l) imply a secular deviation 
given by 

dAB ^ . (37) 

4 Concluding remarks 

We have studied how a small extended body at rest interacts with an incom- 
ing single plane gravitational wave. The body is spinning and also endowed 
with a quadrupolar structure, so that due to the latter property it can be 
thus considered as a good model for a gravitational wave antenna. We have 
discussed the motion of such an extended body by assuming that it can be 
described according to Dixon's model and that the gravitational field of the 
wave is weak, i.e. the "reaction" (induced motion) of a "gravitational wave 
antenna" (the extended body) to the passage of the wave. 

We have found that in general, even if initially absent, the body acquires a 
dipolar moment induced by the quadrupole tensor, a property never pointed 
out before in the literature. As a byproduct, the derivation of the "world 
function" in the weak field limit of a gravitational wave background (absent 
in the literature to our knowledge) has been accomplished and widely used 
throughout the paper. In particular, we have applied it to estimate the or- 
bital deviations between geodesies and the world lines of structured bodies. 
Such deviations could be measured by gravitational wave detectors. Further- 
more, due to the existence of cumulative effects, measuring the deviation of 
the center of mass of the body as induced by the passage of a gravitational 
wave should be favorite with respect to measuring the amplitude of the grav- 
itational wave itself, whose estimate can also be achieved indirectly, since it 
is simply related to the magnitude of the deviation. 

A Limiting cases 

We will examine separately below the complementary cases of an initially non- 
spinning body with a given quadrupole structure as well as that of a purely spinning 
body with a negligible quadrupole structure. 

A.l Body initially not spinning 

In the special case in whi ch the body is initially not spinning (i.e. (tJ = 0), the 
general solutions (|27p -(|29 p reduce to 

ih = ——^ [/()^+ smujTu + 2ql^ Ax {cosujtu — 1)] , 

-1 m 23 a 
J^p = Ax , 



10 



-2 ^ 13 A 
-3 1 12 < 



^foA+ sincuTu + -qo^Ax {cosujtu — 1) , 

-2 1 < 12 . 1 ^ 13 / IN 

!^ = — 2 +'?o smojTt/ — -Ax^o (cosa;T(7 — 1) , 

-3 ^ \ Vi ■ 1 < 12/ -IN 

V = 2"^+'3'o smtJT(7 - -Ax^o (cosa;r(7 - 1) , 

5-^ = ^/o^x sintjrc/ + A+(cos cjTfj - 1) , 

cr^ = -i^+go^(cosa;r!7 - 1) + ^-4xgo^ sintJTLr , 

o-^ = -i^+go^(costjri7 - 1) + ^-4xgo^ sint^ru , (38) 

and the spatial orbit becomes 

ujx = -i/o^+(costjr(7 ~ 1) + igo^Ax (sina;r(7 - ujtjj) , 

= ^^+<7o^(cosa;Tc/ - 1) - iAxgo^(sina;r(7 - ujtu) , 

iLtz = -~^A+ql^ {cos LOTu - 1) - ^Axqo^isinujTu — cotu) ■ (39) 

From Eqs. (|38p it is evident that even initially non-spinning the body starts to spin. 
If the body under consideration is such that q^^ is diagonal we find 

mo . 

m = — — JoA+smuiTu , 

-1 ^ r A ■ ^ -2 n -3 n 

4 mo 

= — foA+ sintjry = — , = Q , — Q , 

4' mo 

= ifoAx sinwrc/ = -24^— , = , = , (40) 
2 y4+ mo 



and the spatial orbit becomes 



ujx — — — /oA+(cos tJTi/ — 1), Loy = Q , cu z = . (41) 



A. 2 Purely spinning particle 

In the case of a purely spinning particle, i.e. with vanishing components of the 
quadrupole tensor, we have: 

m = , 

-1 < 2 3 

Vp = [ctqAx sinujTu + (Joyl+ (cos Ci^Tu — 1)] + ^Ax^Qcro , 
vl = 1: [(^oAx sia-LOTu - cro^+(costi;r(7 - 1)] + ^Ax^lol , 
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= i[((7o)^ - ((To)^]A+ sincjTc/ + o-qo-oAx (cosojrc/ - 1) , 

2 X r 12 2 ~1 

[^xcr(Vo + A+o-q] (cosa;T(7 - 1) , 

l'^ = ^ [A + O-qO-o + AxCTq] shlLUTU 

[Axo-Qcro + A+o-q] (cosa;Tc/ - 1) , 

5^ = , 

o"^ = — yl x ctq (cos a;r(7 — 1) — —A^ag sinujTu , 

o""^ = — yl X ctq (cos cjTy — 1) — — A+ctq sincjT(7 , (42) 



and the spatial orbit 



Ljx = -^[(ctq)^ — {ao)^]A+ {cos LUTu — 1) + o-flcroAx (sintjr(7 — ujtu) , 

2 3 

= -^{A+al + Ax){cosLOTu - 1) - -;f(Axcro + A+){smujru - lotu) 



2 ' ^ " - / 2 



"2 



3 2 

tJ2 = -^(A+ctq + Ax) (cos ojTc/ - 1) - -^{Ay^al + A+){amujru - Loru), (43) 



in agreement with the results of [8]. 

In the special case where the initial spin of the particle is aligned with a trans- 
verse spatial direction, e.g. erg — — a^, the above relations have the following 
limit: 



m = 



1 ' ■ 1 ' 

Up ^ , 1/1 = --AxCTq sinojTc/ , Dp = --A+al{cosLJTu - 1) , 

= ^A+{ao)'^ sin LUTu , v'^ = — ^ Ax ctq sin ajrjj , D^ = —^A+aQ{cosLOTu — 1) , 

cr""^ = , = — iA+(To sintJTLr , — i A x ctq (cos ajT(7 — 1) , (44) 



and the spatial orbit 



LJX ^ -iA+(ao)^(costjr(7 - 1) , 

1 2 

Ujy ^ -Ax(Jo{cOSLJTU - 1) , 

UJ z = — — A+ao(sinajri7 ~ lutu) ■ (45) 
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